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The importance of models with an exact solution for the study of materials with non-trivial 
topological properties has been extensively demonstrated. Among these, the Kitaev model of a 
one-dimensional p-wave superconductor plays a guiding role in the search for Majorana modes in 
condensed matter systems. Also, the sp chain, with an anti-symmetric mixing among the s and p 
bands provides a paradigmatic example of a topological insulator with well understood properties. 

There is an intimate relation between these two models and in particular their topological quantum 
phase transitions share the same universality class. Here we consider a two-band sp model of spinless 
fermions with an attractive (inter-band) interaction. Both the interaction and hybridization between 
the s and p fermions are anti-symmetric. The zero temperature phase diagram of the model presents 
a variety of phases including a Weyl superconductor, topological insulator and trivial phases. The 
quantum phase transitions between these phases can be either continuous or discontinuous. We 
show that the transition from the topological superconducting phase to the trivial one has critical 
exponents different from those of an equivalent transition in Kitaev’s model. 


I. INTRODUCTION 

Since the first strong experimental evidence of 
Majorana fermions [1] in a hybrid superconductor- 
semiconductor one dimensional system, the search for ex¬ 
otic states supporting Majorana fermions has attracted 
increasing interest in condensed matter physics. Recent 
observations have reinforced the existence of Majorana, 
specially the one made in ferromagnetic atomic chains on 
a superconductor (SC) J2|. Anomalous behaviour on this 
experimental evidence |3) indicates that the appearance 
of Majorana may have yet unknown sources. Also some 
technical difficulties such as highly localized states com¬ 
pared to the material parameters |4] or high temperatures 
that prevent to have access inside the gap [5] have left 
the existence of Majorana inconclusive. The running for 
the experimental discovery of Majorana is well described 
in 0. 

It is well known that the Kitaev model [Ml - anti¬ 
symmetric pairs of spinless fermions in ID - exhibits a 
non-trivial topological phase with Majorana modes at the 
ends of a p -wave superconducting chain. The excitations 
at the ends of the chain depend on the quantum state 
of the system, which in turn is determined by the ratio 
/i/2£, between the chemical potential (i and the hopping 
t. If \n\/2t < 1 the chain is superconducting with non¬ 
trivial topological properties. This weak pairing phase 
presents Majorana fermions at its ends. Otherwise, if 
|/i|/2 1 > 1 the chain is in a strong coupling supercon¬ 
ducting phase with trivial topological properties and has 
no end states 0. 

The importance of the mixing of sp bands for topo¬ 
logical insulators has already been pointed out in differ- 
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ent contexts, including that of the spin quantum Hall 
effect nni and a cold atom version of the sp-chain, the 
sp- ladder[11]. The sp mixing is in a special class that 
mixes orbitals with angular momenta that differ by an 
odd number. This implies the anti-symmetric prop¬ 
erty V (—r) = —V (r) or in momentum space V(— k) = 

— v(k) m ■ 

In addition, it was recently shown Q3HIS] an inti¬ 
mate relation between a two band insulator with anti¬ 
symmetric hybridization and the Kitaev model, as re¬ 
gards to the topological properties and their end states. 
By tuning the parameters of the ID sp chain, the sys¬ 
tem can be driven, through a topological quantum phase 
transition, from a trivial to a topological insulator. As a 
result they found two Majorana zero modes at the ends 
of an insulating chain. In the search for Majorana modes, 
Kitaev’s model is a most clear example of the importance 
of exactly soluble models as guides in this difficult path. 

Here we consider a different model of a p-wave super¬ 
conducting chain that also can be solved exactly. This 
study will throw further light on the role of topology 
and different type of interactions and symmetry break¬ 
ing terms as conditions for the existence of Majorana in 
superconducting wires. We consider a two-band model 
of spinless fermions in a chain with inter-band attractive 
interactions and an anti-symmetric hybridization. The 
model is exactly soluble and we obtain its zero temper¬ 
ature phase diagram. For a fixed small value of the 
mixing, as the ratio /i/2 1 of the model increases, such 
that |/i|/2 1 > 1 , there is a topological quantum phase 
transition from a gapless (topologically non-trivial) to a 
gapped (trivial) superconducting phase. The phase dia¬ 
gram resembles that of Kitaev’s model, however the na¬ 
ture of the topological phases, as well as, the topologi¬ 
cal transition are distinct. In our model the non-trivial 
superconducting phase has gapless Fermi points. These 
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gapless points have the characteristic of Weyl fermions 
in 3D systems mm, as they have a non-degenerated 
linear dispersion relation and appear or disappear in 
pairs, only when two Fermi points unite. The conserva¬ 
tion of the topologicalcharge associated with these Fermi 
points m confers a non-trivial topological character for 
this phase. Furthermore, we show that the topological 
quantum phase transition at (\fi\/2t) c = 1 is in a differ¬ 
ent universality class of that of Kitaev’s model. Next, 
fixing the chemical potential, say at /i = 0, and increas¬ 
ing the hybridization the system is driven from the Weyl 
superconducting phase (WSC) to a topological insulator 
through a first order quantum phase transition. Topolog¬ 
ical phase transitions are known to produce anomalies in 
thermodynamic quantities[20] and we obtain these here, 
with special emphasis on the behaviour of the compress¬ 
ibility. 

In section[TI]we define the model and analyse the super¬ 
conductivity stability. In section m we show the phase 
diagram for the superconducting stability. We also dis¬ 
cuss the nature of the transitions and show the localiza¬ 
tion of end states in a finite system. In section [IV] we 
calculate topological invariants of different phases, and 
discuss the topological properties of the model in Majo- 
rana basis. Finally, in section [V] we present some conclu¬ 
sions and highlight the main results. 


II. DEFINING THE MODEL 

We consider a two-band problem with hybridisation 
and triplet inter-band superconductivity in ID, i.e., a 
chain with two orbitals per site, with angular momenta 
differing by an odd number, let’s say p and s. The pairing 
between fermions on different bands (inter-band) is al¬ 
ways p-wave kind, in the sense that the pairing of spinless 
fermions is anti-symmetric. The problem can be viewed 
as a generalization of Kitaev’s model to two orbitals and 
only interband pairing. We also have the anti-symmetric 
hybridisation term that, under some conditions, was 
shown to be responsible for topological phases [13415]. 
The simplest Hamiltonian in the momentum space that 
describes those types of superconductivity and hybridiza¬ 
tion can be written as 

U = ( c\c k + f4p fe ) + 2 1 cos (k) (p{p k - 4c fc ) 

k 

-iA sin ( k ) c\p^_ k + iV sm(k)c\p k + h.c.j , 

( 1 ) 

where fi is the chemical potential, A is the sp pairing am¬ 
plitude, and V is the anti-symmetric hybridization am¬ 
plitude. Note that the hopping amplitude t has different 
sign in each band, representing particles for the orbital 
s and holes for the orbital p. We can write the same 
Hamiltonian using the Bogoliubov-de Gennes (BdG) rep¬ 
resentation as 


H = Y, C l'HkC k , ( 2 ) 

k 

with C\ = ( c\p\c- k p-k^ and 

Bk = /il\ 0 zz T ^k^yx VkJ^zy> (3) 

where T ab = r a (g) r b , V a, 6 = x,y,z, and r x , y:Z /r X:y , z are 
the Pauli matrices acting on particle-hole/orbitals space, 
respectively, and ro = To are the 2x2 identity matrix. We 
have defined Sk = 2£cos(£;), V& = Vsm(k), and A/~ = 
A sin ( k ). 


A. Energy spectrum 


Since a topological phase transition only occurs when a 
gap closes, looking for gapless points on the energy spec¬ 
trum may indicate this transition. The model considered 
here has the following energy dispersion relations, 


E(k ) = ±\J Zi ± 2v^2, (4) 

where Zi = A(k) + B (k) and Z 2 = A(k)B(k), with 
A(k) = + and B (k) = A 2 +/i 2 . Looking for gapless 

points (E(k) = 0) the possible solution is A(k) = B (fc), 
i.e., 


M 2 = 4 + n 2 -A l 

= (2 tf cos 2 (k) + (V 2 - A 2 ) sin 2 (k) 


(5) 


We will analyze the equation above more deeply in sec¬ 
tion m First we would like to highlight the case with 
no hybridization, V = 0, in which the system is always 
gapless whenever \fi\ < 2 1. The existence of these gapless 
modes represents a substantial difference between this 
and the Kitaev model, see figure [1] We will see in the 
next section that even in this non-gapped region the sys¬ 
tem shows superconductivity. On the other side, when 
\fi\ > 2£, the system is fully gapped but superconductiv¬ 
ity is still present up to \fi\ <4 1. 

Deep inside the gapless phase the crossings between 
bands have a linear dispersion relation @0 and define 
Dirac nodes. Furthermore, we note that the bands are 
non-degenerate and the nodes appear and disappear only 
when two nodes are combined, as one can see com¬ 
paring figures & and [TJd. This is a characteristic of 
Weyl fermions in 3D or 2D SC[T8] and in topological 
superfluidity [16, 17]. In this sense, the model here pre¬ 
sented can be called ID Weyl SC. 
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Figure 1. Dispersion relations for a pure sp superconducting system, without hybridization, for A = 1.5 and V = 0. The 
system supports gapless excitations when (\g\ /2t) < 1 and is gapped otherwise, specifically, we set (|/i| /2t) = 0.25 (left) and 
(\fi\ /2t) — 1.5 (right). In this work we show that both regions have self-consistent solutions for the superconductivity. The 
central panel shows the zero modes annihilation; note that it happens in pairs of nodes. 


B. Self-consistent equations for the 
superconductivity and the occupation number 


We may calculate the self-consistent inter-band super¬ 
conducting order parameter, A, from the gap equation 

A = -j- ^isin(fc) (p_ ft Cfc) , (6) 

k 

where, g is the attractive energy between the spinless 
fermions, L is the length of the chain, and the correla¬ 
tion function ( p- k Ck) is obtained from the fluctuation- 
dissipation theorem[2T] (a similar calculation was re¬ 
cently done in [13): such that 


(p~ k c k ) = 2tt / i(( c k’P-k)) r - ((c k ,P-k)) a ]dui; 

(7) 

the Green’s functions ((c/ c ,_p_/ c )) r,a (retarded and ad¬ 
vanced) are obtained from the Greenian operator [221124] , 
i.e., (( G)) k = (cjI 4 x 4 — Hk)~ , and f (cj) is the Fermi 
distribution. If we proceed with the calculations, see ap¬ 
pendix [Aj the gap equation at T=0 becomes 


where 


1 1 4 sin 2 (k) 

g L Z-^ (uq + oq) 



if B(k)> A (fc), 
otherwise. 


( 8 ) 


( 9 ) 


Also, uq and are the eigenvalues of the Hamiltonian, 
such that, 


LUi(k) = E (fe) + — \JZi + 2 y/Z^, 

&2{k) = E (k)_ = \Jz i - 2y/ Zv 2 . 

We can verify the stability of the superconduct¬ 
ing phase calculating the parameters A and g self- 
consistently from the gap and the occupation number 


equation given by, 


n = n s + n p = [( 4 C *) + (PkPk 


( 11 ) 


with n s and n p the occupation numbers for the 5 and p 
bands, respectively. For the model considered here, the 
above equation is 



a 

(uji + CJ 2 ) 


4 + 


1 

2 


( 12 ) 


where 0 < n < 2 is the total occupation number per site 
of the chain. 


III. PHASE DIAGRAM 


The solution of the coupled self-consistent equations 
for the gap and the chemical potential is complicated 
by the constraints of the sums in momentum space (Eq. 
(|9|). In this section we present the phase diagram of our 
model system obtained directly from a numerical solution 
of the BdG equations fixing the chemical potential. The 
Hamiltonian defined in Eq. <© can be solved using BdG 
transformations as 


C k = Y \ U n,kln,k + A 


Pk = 




[<fc7n,fc + (<fc) 


y 

in, — k 


(13) 


This transformation diagonalizes the Hamiltonian in the 
form 


n k 


<,fc \ 
< k 

\<-j . 


/ u: 


= E n 


n,k \ 
P 

X n,k 


V n,-k 

\ V n,-kJ 


(14) 


where E n are the energy eigenvalues and the wave func¬ 
tions spinor are the eigenstates. The self-consistent solu¬ 
tion implies that the pairing can be obtained using 

A = sin(fc) (( c- k p k ) + (p-feCfe)). (15) 

k 
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At zero temperature, using the representation of the A. Nature of the transitions 

fermionic operators in terms of the wave functions and 

the Bogoliubov coefficients we may write g WSC-TI Half-filling 


A 


2 ^ 

L 


k n 




(<-*)'+<*<* 



In Fig. we show the numerical results for the or¬ 
der parameter A as a function of the chemical potential 
and hybridization for a fixed value of the attractive in¬ 
teraction g = 1.7. All quantities are normalized by the 
hopping term t. In Fig. |2j: we show the gap for exci¬ 
tations for the same range of parameters. The results 
in these figures allow us to obtain the zero temperature 
phase diagram of the system shown in Fig. It- 

In agreement with our previous discussions we find 
a gapless superconducting phase for (\g\/2t) < 1 and 
V/21 < V c (g)/2t = [(g/2t) 2 + (A 0 /2t) 2 ] 1 / 2 named Weyl 
SC (WSC) in the phase diagram. We note that at the 
transition (V = V c (g)) the gapless points always occur at 
fc = ±7r/2. The quantity Ao = Ao(/i) is the value of the 
order parameter at V = 0 for a given chemical potential 
value g. 

For (\g\ /21) > 1 and V < V c (g) the system presents 
a gapped superconducting phase with trivial topological 
properties similar to the strong coupling superconduct¬ 
ing phase of Kitaev’s model. In this phase, named SC 
in the phase diagram, the order parameter vanishes con¬ 
tinuously as the chemical potential increases. For a fixed 
V < V c (g), the range of this phase for increasing g/21 
depends on the strength of the attractive interaction g. 

On the other hand for (\g\ /21) < 1, but for V > V c (g) y 
there is a gapped non-superconducting phase, that cor¬ 
responds to a topological insulator (TI), as will be dis¬ 
cussed below. This phase is characterized by zero energy 
modes localized at the ends of the chain for g = 0. There 
are also localized modes if g ^ 0 that have finite subgap 
energy. 

Notice that the conditions for the existence of a gap are 
given by Eq. <§• For instance, in the case of strong hy¬ 
bridization and weak or no superconductivity, such that 
l + (A/2£) 2 < (V/2t) 2 , the system becomes gapless when¬ 
ever (g/2t) > 1 and V > V c (g). This corresponds to the 
phase M in the phase diagram of figure §> which is a 
normal metallic or insulating phase (not shown in the 
figure) depending on the occupation number. 

In order to clarify the understanding of the phases dis¬ 
cussed above, in figure [3] we plot the energy spectrum in 
different regions of the phase diagram of figure^. These 
figures illustrate the cases of appearance of gapped or 
gapless superconductivity. 


Let us first consider the case of half-filling bands 
(n = 1), where the chemical potential is fixed at g m 
0. The constraint in equation ^ now reads A 2 > 
(e^ + V^). Important points in momentum space cor¬ 
respond to those wavevectors where this inequality be¬ 
comes an equality, i.e., 


tan 2 (fco) = A 2 2 ^ F 2 - (17) 

At these points the system becomes gapless and they 
characterize the Weyl points. From this result, we can 
immediately see that there are no gapless nodes when 
V > A. It is easy to see also that with increasing hy¬ 
bridization the Weyl points collapse at fco = vr/2 = kp 
for V = V c = A at a discontinuous quantum phase tran¬ 
sition from the WSC to the TI phase where the super¬ 
conducting order parameter drops to zero. This collapse 
of superconductivity is associated with the appearance 
of zero energy modes exactly at the Fermi surface kp, of 
the half-filled system. In the superconducting side before 
the transition, the order parameter (when g = 0) is given 
by the gap equation, 


1 

9 


4 

\A[L 


k 0 

J2 sin ( fc ) • 

k =0 


(18) 


where fco is the largest momentum value that contributes 
to the superconductivity. Notice that this equation has 
no trivial analytic solution since fco depends on A. At 
the transition, for V c = A the momenta fco = ±7r/2 and 
A 0 = 4g/27r before dropping abruptly to zero at the TI 
phase. 


2. WSC-SC 

We now investigate the transition from the non-trivial 
topological superconductor to the trivial one by increas¬ 
ing the chemical potential at fixed hybridization. Let us 
for simplicity consider the case of V = 0. The WSC-SC 
transition occurs for ( g/2t) c = 1 as shown in figure [2]:. It 
is associated, as can be easily checked with the collapse 
of two Weyl points at the center of the Brillouin zone 
(fc = 0) and at its extremities (fc = ±7r). Expanding the 
dispersion relation of the excitations close to fc = 0 and 
(>u/2t) c = 1, we get, 


^(fc) = 2t^(l-|)2 + (^)2fc4. (19) 

We have omitted the fc 2 term, since its coefficient is pro¬ 
portional to (1 — g/2t) and vanishes at the quantum 
topological phase transition at ( g/2t) c = 1. Then, at 
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Figure 2. The superconducting amplitude A (left panel) and spectrum gap (right panel) for different values of V and g, with 
g — 1.7. The middle panel shows all the different phases, that are: trivial superconducting gapped phase (SC), metallic (M), 
gapless Weyl superconductor (WSC) and topological insulator (TI). 




Figure 3. Energy spectrum in momentum space for different points on the phase diagram shown in figure [2] The up left 
panel shows a system with hybridization and superconductivity when V < A, and the chemical potential is fine tuned to 
zero (specifically, t = l, /z = 0, V = 0.5 and A = 1.5); the up right panel illustrate the case where there is no self-consistent 
superconductivity (specifically, t — 1, g = A = 0 and V = 1.5). The down panels show the cases where there is superconductivity 
but the chemical potential is driven out from zero. The left plot has the same conditions that the first case but with g = 0.5, 
and it illustrates the case with gapless superconductivity. On the other hand, the right panel was set to same conditions as 
before but g = 3 and A = 1, and it illustrates the full gapped superconducting region. 


the quantum critical point, the spectrum of excitations 
uj2{k)(x k 2 , which allows to identify the dynamical ex¬ 
ponent z = 2 for this transition. On the other hand 
at k = 0, the gap uj 2{k = 0) = (g/2t) c — (g/2t), van¬ 
ishes linearly at the quantum critical point with a gap 
exponent vz = 1. The critical exponents v = 1/2 and 
z = 2 show that the quantum phase transition from 
the topological to the trivial superconducting phase in 
the inter-band model is in a different universality class 
from that of the Kitaev model. In the latter at the 
QCP, ( g/2t) c = 1 , the dispersion is linear implying a 
dynamic exponent z = 1. Since the gap vanishes lin¬ 
early also, we get for the correlation length exponent the 


value v = 1 (see Ref. El)- These different values of the 
critical exponents imply distinct behaviour for the com¬ 
pressibility of the two models at the topological quan¬ 
tum phase transition inside the superconducting phase. 
The compressibility close to this transition is given by, 
Xc = d 2 f/dg 2 oc \((g/2t) c — (g/2t)\~ a where / is the free 
energy density. The exponent a is related to the corre¬ 
lation length and dynamical exponents by the quantum 
hyperscaling relation (25H27] , 2 — a = v{d + z). It can be 
easily verified that while for the intra-band Kitaev model 
a = 0, which is generally associated with a logarithmic 
singularity [28], for the inter-band model a = 1/2 imply¬ 
ing an even stronger singularity for the compressibility 
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at the topological transition. Indeed in our model the 
topological transition is in the universality of the Lifshitz 
transition m- Notice that this is a purely topological 
quantum phase transition, since both phases are charac¬ 
terised by the same order parameter. In spite of this, they 
have singularities described by critical exponents which 
obey the quantum hyperscaling relation [25]. Although 
the usual Landau approach of expanding the free energy 
in terms of order parameters that become small close to 
a continuous phase transition is of no use here, the renor¬ 
malisation group still provides an adequate description of 
this critical phenomenon m- 


B. Fermi velocity 

We may calculate the Fermi velocity at the Fermi 
points fco expanding the energy spectrum in equation 
0 in their vicinity. The spectrum becomes E (fc) = 
E (fco) + vp (fc — fco) + O (fc 2 ), where vp = Vf (fco) is the 
Fermi velocity. For the general condition in equation 0 
the Fermi points are given by, 


fco = arcsin 

where, of course, the term inside the brackets must be 
within the range [0,1]. The Fermi velocities were ob¬ 
tained for three different situations and are shown in 
Fig-EI The up row shows the variation of fco for a sys¬ 
tem without hybridization (V = 0), for a fixed A = 2, 
then for fixed A = V = 2, respectively from left to right. 
The down row shows the Fermi velocities for each case, 
corresponding to the upper plot. 


I (2 ty-iE 

A 2 - 1/2 + (2 tf 


( 20 ) 


C. Energy spectrum in real space 


In order to find the energy spectrum in real space 
through the BdG transformation, we write the Hamil¬ 
tonian in the form, 

H = C^HC, (21) 


where 

c = (c 1 Pi c\ p\ ■■■ c N p N c ] N p ] N ) T (22) 
and the matrix H is defined as 

(Hu ••• Ui N \ 

H : •. : , (23) 

\Hni ••• 'Hnn) 

and is comprised by the following (4 x 4) interaction ma¬ 
trices 


Er,r = fjT z0 , 

^r,r+1 = tY zz T iVY Z y, 

y~Lr,r —1 = tYzz T i^Yyx iVYzy, 

k~L r y =0 V r' / r, r + 1 or r — 1. 


( 24 ) 


The BdG transformation, 

C, = £ [u s ,n( r hn + <„(07n] - 

n 

Pr = [^p,n( r )7n + V p,n( r )lli\ > 

n 


(25) 


diagonalizes the Hamiltonian, T-L = Eq + £ Enlfnlfm 

n 

such that, 


U ] HU = E , (26) 

where U is formed by all the BdG coefficients u s , v s , u p 
and v p , and has the property to be unitary U^U = I. The 
matrix E is diagonal and contains the energy spectrum 
(E n ) of the system. 

We have calculated the energy spectrum for a chain of 
L = 100 sites with two-orbitals per site and inter-band 
interactions in the presence of hybridization. The spec¬ 
trum consists of 4 L energies. We have checked that this 
size is large enough to prevent finite size effects. We set 
the chemical potential to zero fi = 0, and take the hy¬ 
bridization strong enough (V > V c or A = 0), such that 
the system is in the TI phase of the phase diagram in Fig. 
It- We can see in Fig. 0 the appearance of zero energy 
fermionic modes (see below). We also checked that these 
zero energy modes are localized at the ends of the chain, 
see Fig. 0- This gapped insulating phase share the same 
properties of the topological insulating phase found in a 
normal sp chain [IMl]. In this situation the fermionic 
modes resemble the Majorana zero modes, as will be dis¬ 
cussed in the next section. 

Next, we remove the chemical potential from zero and 
keep (V > V c ), but such that the system remains in the 
TI phase. We now find that there are two (plus two 
particle-hole symmetric) energies in the spectrum dis¬ 
placed from zero energy, corresponding also to localized 
edge states but of finite energy and longer spatial extent. 

A more intriguing situation happens when we induce 
superconductivity on this TI phase and other energies 
(for a small range of parameters) displaced from the spec¬ 
trum appear, see figure [5] down left panel. The system 
now have 4 localized states (plus 4 particle-hole symmet¬ 
ric) consisting of two double-degenerate states. We also 
show that all those particular energies are localized in 
the end of the chain, but with lower occupation number 
each, see down right panel of figure [5] 


IV. TOPOLOGICAL INVARIANTS AND EDGE 
MODES 

A. WSC - topological invariant 

The winding number is a proper topological invari¬ 
ant that classifies the topological phase of a gapped ID 
system. In the gapless superconducting phase, it is not 
possible to calculate this by conventional methods since 
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Figure 4. Up row shows the positive Dirac points for a half-filled system (fi = 0) without hybridization (V = 0), for a system 
with fixed superconductivity value (A = 2), and for fixed hybridization and superconductivity (A = V = 2) when increasing 
the chemical potential, respectively from left to right. The down row shows the Fermi velocity for each case, corresponding to 
its upper plot. 


Energy 



«o(i) 



Energy 



n o0) 
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Figure 5. Energy spectrum in real space for a chain of L — 100 sites. In the first row, the left panel shows a hybridised system 
with no, or weak induced, superconductivity when the chemical potential is fine-tuned to zero, n = 0. The isolated point is 
four-fold degenerated. On the right, the panel shows the localisation of the Majorana zero modes. In the second row, the left 
panel shows the same system, with induced superconductivity, but the chemical potential has its value near to zero. Those 
isolated points are two-fold degenerated. The right panel shows the occupation number on each lattice site for one of those 
isolated energies. It is clear that they are localised states at the edge. 


there are zero energy points that cannot be avoided in 
one dimensional systems, or the sum over the Brillouin 
zone gives a vanishing winding number since the Fermi 
points appear in pairs and their contributions cancel out 


[29] , On the other hand, let’s look closer to one of the 
linear dispersion relations (see fig. [3| that crosses the 
zero energy at some point k = k$. In this region the 
Hamiltonian with Weyl nodes in ID can be reduced to 


























describe the two Bogoliubov bands that cross zero en¬ 
ergy. The reduced low energy part of the Hamiltonian 
may be expanded in terms of Pauli matrix such that 

H{k)= ]T di(k)(Ji. (27) 

*= 1 , 2,3 

Imposing inversion and charge conjugation operations 
m we end up with only one independent term, for ex¬ 
ample, d z (k) ~ (fc — fco), where ±fco are the two gapless 
points and the Hamiltonian is simply 

H = (k - fc 0 ) (J z 

= <W (28) 

For a positive chemical potential the energy at the Fermi 
level is positive, which yields the eigenstate of matrix o z 
to be +1 or —1 depending on the momentum fc, such that 
a ^ if k < fco and <j z i/j = otherwise. In this 
frame the winding number can be calculated as [30] 

W = sgn (4) (ip~ 1 a z ip) . (29) 

Therefore, W = +1 indicates a non-trivial phase with 
topological excitations, the Weyl nodes. Note that if the 
system is not Weyl-like the Hamiltonian cannot be writ¬ 
ten in terms of one independent component and the wind¬ 
ing number value is no longer preserved. Furthermore, 
the above result is valid not only in the Pauli basis but 
also for higher order matrices, provided that codimension 
is zero [29], such as the 4x4 Hamiltonian here considered. 

Let us look for a special case of half-filling bands and 
no hybridization, fi = V = 0. In this situation we are 
in the WSC phase where the Weyl fermions appear at 
the momenta k = ± 7 r /4 and fc = =L 37 r/ 4 , when A « 2 1. 
If we make a basis rotation on eq. (p| such that C\ — 
[c\p-kC-kV^ it is easy to see that it can be decoupled 

in two 2 x 2 Hamiltonians. Near to the Weyl point fco, 
one of these two Hamiltonians, e.g., the one for the basis 

( 4 p- fe ), is 

H = (a 0 (fc - fc 0 ) + mo) cr 0 + (a(fc - fco) + m) a y , (30) 

with <ao = mo = 2t/ \/( 2 ) and a = m = A/y/(2). Disre¬ 
garding the mass term mo which leads to a shift on the 
energy, we have only one independent term, the mass 
term m, which cannot produce a gap by itself on the 
spectrum. 


B. TI - Majorana modes 

In order to clarify the existence of Majorana modes 
in our model, we write the Hamiltonian, Eq. [l] in real 
space. This is given by, 

u = E {~P ( C k +PiPi) + 1 {PiPi+l - 4 c *+i) 

+V (c\p i+1 - 4 +1 pij + A (pjCi +1 - Pi+iCi) 

+h.c.}. (31) 


This can be written in terms of Majorana operators, a^ r , 
Pa,v and via the relations, 


C r — ^ (^ B,r T ^^A,r) and p r — ^ (@B,r T ' 

( 32 ) 

Now, we perform a second transformation on Majorana 
fermions - we call them unconventional hybridized Ma¬ 
jorana fermions [311 - such that 


&A,r 


@A,r — 


jXr + 7 A,r 

2y/(V-A)’ 


&B,r 


7 B,r + 7 B,r 

2 \/ W + A) 


(33) 


7tp - 7 A . 


2 y (v + a) 

the result is the following, 


and pBr = 


7 R, 


' 7 ba 


2VW -A)’ 


U ' = i E {(! - C ‘) fer'rtr+l ~ ^A,rlB,r+l) 
r 

- (1 + Ct) (lB, r lA,r+l ~ 7A,r7B,r+l) 

(7B,r7t,r + 7s, r 7^,r) } > (34) 

where C t = tj\/ ( V 2 — A 2 ) and = \ij\] ( V 2 — A 2 ). If 
we go to the limit p = 0 and take C t = 1 , such that 
the system is in the TI phase in Fig. (V > A), the 
Hamiltonian reads, 

. N—l 

^ = — 2 {^B,r^A,r+l ~ ^A,r^B > (35) 

r=l 

which couples Majorana fermions only at adjacent lattice 
sites. Proceeding with the same analysis as 0, we may 
easily see that the Majorana modes 71 = 7^4 l5 72 = 7 ^ 1 , 

73 = 7 b tv’ anc ^ 74 = 7 a at are n °l present in the above 
Hamiltonian, it means that they have no cost of energy 
to be added to the system; they are called Majorana zero 
modes. In the present case we have two Majorana zero 
modes on each end, and they can be combined to form 
one ordinary fermionic operator at each end as 

h = \ (7i + * 72 ) and f N = | (73 + * 74 ), (36) 

or, naively, we may think of two highly non-local 
fermionic operators such as 

/ = \ ( 7 i + * 73 ) and f = \ (72 + «74) • (37) 

One may note that this result agrees with the conclusion 
on figure [5] up right panel. Moreover the Majorana zero 
modes persists out of the fine-tuned C± = 1 , or t 2 = 


that there is no gap closing for this range of parameters. 

When /i 7 ^ 0, the topological character is preserved, in 
the sense that we still have localized states at the ends of 


V — A' , provided /i = 0, since we know from section IIA 
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the chain, but the Majorana zero modes are not robust 
such that they acquire a finite energy, see down left panel 
of figure [5] In this situation, or when general referring, 
we call them by fermionic modes, instead of Majorana 
ones. 


C. TI - topological invariant 

The non-trivial topological character of the TI phase, 
can be shown by calculating the winding number for the 
special case /x = 0. In this region of the phase diagram we 
have A = 0 and the 4x4 Hamiltonian can be decoupled 
in two 2x2 Hamiltonians, such as 

'H = —Vk<J y — Sk&z- (38) 

This equation can be rewritten as the Hamiltonian of a 
spin 1/2 in a ^-dependent magnetic field, 

H = —h (k) • cr, (39) 

where h (k) = ( h x , h yi h z ) = (0, — 14 , —£k) with the prop¬ 
erties h Xr y(k) = —h x , y (—k) and h z {k ) = h z (—k). The 
winding number v is obtained as the product of the signs 
of the magnetic field on the center and at the extreme of 
the Brillouin zone, i.e., 

v = sgn(h (k = 0))sgn(h(fc = i r)). (40) 

Since h (k = 0) = (0, 0, 2 1) and h(fc = tt) = (0, 0, —2£), we 
get v = — 1, which characterizes the non-trivial topolog¬ 
ical character of the TI phase along the line /x = 0. The 
topological nature of this phase is associated with the 
existence of zero energy modes at the ends of the chain, 
as discussed above in sec. IIVBI and also calculated and 
shown in figure [5] 

As /x increases, we observe from the numerical solution 
that the end modes persist on the chain but acquire a 
finite energy. As /x/2£ = 1, where the gap of the TI phase 
vanishes, they merge with the continuum of excitations. 

If one calculate the winding number (by usual meth¬ 
ods) for the whole system, the 4x4 Hamiltonian including 
A, it shows itself trivial. The topology is hidden by the 
charge conjugation (or particle-hole) symmetry imposed 
to the system. Some attempts to calculate the winding 
number using new methods were proposed to uncover this 
kind of topology. |32| . On the other hand, the topologi¬ 
cal character of the whole TI phase is guaranteed since it 
is adiabatically connected with the topological case just 
shown (when A = fi = 0). 

V. CONCLUSIONS 

We have studied in this work a ID sp-chain with at¬ 
tractive inter-band interactions and anti-symmetric hy¬ 


bridization due to the different parities of the s and p or¬ 
bitals. The latter was shown to be responsible for the ap¬ 
pearance of topological phases [musij in non-interacting 
sp-chains. We have shown that this model presents a rich 
phase diagram including non-trivial topological phases. 
It is interesting to compare it with the Kitaev model 
which also has an exact solution. In both models there is 
a weak coupling superconducting phase with non-trivial 
topological properties. However, while in Kitaev’s model 
this phase is gapped, in our model it has Fermi points 
with gapless excitations. 

We have studied the quantum topological phase tran¬ 
sition between the weak coupling, non-trivial to the triv¬ 
ial, strong coupling superconductor and found that this 
transition in our model is in a different universality class 
from that of Kitaev’s model. In the strong coupling limit, 
the superconductivity disappears if fi is very large. We 
have also shown the existence of a discontinuous quan¬ 
tum phase transition from a Weyl superconductor to a 
topological insulator with increasing hybridization. This 
is caused by the appearance of a zero energy mode ex¬ 
actly at the Fermi surface of the normal, non-interacting 
system. 

We have shown that in the phase diagram of the 
present model there is a topological insulating phase, 
with zero energy fermionic modes at the ends of the 
chain. This phase has been characterized by calculat¬ 
ing its winding number and the zero energy modes have 
been found both analytically and numerically. 

The importance of models with exact solutions in the 
theory of topological matter has been now widely recog¬ 
nized. Besides throwing light in many exotic properties 
of these materials, they serve as guides for obtaining new 
types of excitations which are protected by topology. The 
present model, which as we have shown exhibits a rich 
variety of phases and different types of phase transitions, 
has many new features that allows for a deeper under¬ 
standing of topological systems. 
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Appendix A: Gap equation with hybridization 

In order to demonstrate the result in Eq. <© we start 
calculating the gap equation 

A = -j- ^zsin(fc) (p-fcCfe), (Al) 

k 

that we solved using the fluctuation-dissipation theorem 

(p-k c k) = ^ J /(w) [{{c k ,p- k )) r - {{c k ,p- k )) a ](Lj, 

(A2) 

where / (u) is the Fermi distribution. The retarded and 
advanced Green functions are obtained from the Gree- 
nian ((G)) k = (?G,) _ \ with TL U cjI 4X 4 - 74, such 
that 


m) k = 


Fk (&) 

det (Hu)' 


(A3) 


With the same basis used in equation ©> the Green func¬ 
tion ((ck,P-k}) is the (4,1) element of the Greenian ma¬ 
trix, which can be written as 


(F k — i^kF k (m) 


= — iAk 


(uj ^ e fc ) 2 - /i 2 - + V k (A4) 


If one put this into the Eq. (A2), proceed with the in¬ 


tegral calculation and takes the zero temperature limit, 
after some calculations, will find that mean value of the 
operators is 


(P-kCk) 

where 


— f F k (— UJ2) _ Fk (—^ 1 ) 1 
2 (<jj\ - u%) \ cj 2 cui j ’ 

(AS) 


(k) = \fzf+ 2 yfz ~2 and w 2 (fc) = \fzf~ 2vTsi 

(A6) 

with 

^ = el + V k 2 +At+p 2 (A7) 

and 

Z 2 = (Ai + S)(V> + el). (A8) 
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If we put this result into the gap equation we find that 

1 _ 1 sin 2 (k) f {—(^ 2 ) F^ (— u)\) 1 

4 g L “ 2 (wf — wf) \ w 2 wi J ' 

(A9) 

We can rewrite the right hand side of the equation above 
as 


1 f Fk (-^2) _ Fk (~^i) \ _ 

W 2 LOi ) 

= J^TW) {ssn[A{k) ~ B{k)] ~ 1} ' 


(A10) 


We may write this result in the compact form 


1 f Fk (—^ 2 ) 
(wf - wf) \ UJ 2 


Fk (-wi) 


LOi 


-24 

(wi + 0J 2 ) ’ 

(A12) 


where 



if B(k)>A(k), 
otherwise. 


(A13) 


Such condition shows that the gap equation to this model 
is 


where 

A ( k ) = (V 2 + 4) and B (, k) = (A 2 k + ^ 2 ). 

(All) 


1 1 ^ ^ sin 2 ( k ) 

g L F (loi + uj 2 ) k 


(A14) 



